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A microscopic mean-field description of lieavy ions fusion is performed in the framework of 
the Time-Dependent Hartree-Fock (TDHF) theory using a Skyrme interaction with the SLy4d 
parametrization. A good agreement with experiments is obtained on the position of the fusion 
barriers for various total masses, mass asymmetries and deformations. The excitation function of 
the ^''0-|-^''*Pb is overestimated by about 16% above the barrier. The restriction to an independent 
particles state in the mean-field dynamics prevents the description of sub-barrier fusion. Effect of 
t^^ ■ transfer on fusion is discussed. 

o: 

CN . I. INTRODUCTION 

>: 

_ O ■ Description of nuclear reactions is very challenging, especially at energies around the fusion barrier generated by 

| /' , ' the competition between Coulomb and nuclear interactions. It has been established that the structure of the collision 
partners may affect strongly the reaction mechanisms in this energy domain, as, for instance, the fusion cross-sections 
(for a review, see e.g. Ref. [1]). 

^-1, Upcoming exotic beams facilities at energies of a few MeV/u like SPIRAL2 will allow studies of the interplay 

2 ■ between reaction mechanisms and exotic nuclear structures such as haloes, neutron skins, high isospins... It is 
I , therefore recommended to treat both structure and dynamics within the same formalism. This is the case of fully 
Q ■ microscopic approaches such as the Time-Dependent Hartree-Fock (TDHF) theory proposed by P. A. M. Dirac in 
^ 1930 [2]. This time dependent version of the well known Hartree-Fock (HF) theory gives a self-consistent mean-field 

Ph ■ description of nuclear dynamics [3, 4]. The success of the first HF calculations based on the Skyrme interaction 
[5, 6] led to tremendous activities to describe nuclear structure within mean-field based approaches (see Ref. [7] for a 
review) . 

J> ' Early TDHF calculations has been done with the seek for a description of the dynamics of nuclei as good as their 

^^ , static properties [8-16]. They used various symmetries and simplified Skyrme interactions to reduce the computational 
Cn . time. The increase of computational power allowed realistic TDHF calculations of nuclear collisions in 3 dimensions 
with full Skyrme interactions in the last ten years [17-23]. 

In this paper, we present a TDHF study of nuclear fusion. In a first part, we recall the formalism and detail the 
calculation. In a second part, we calculate fusion barriers for several systems with various total masses and mass 
asymmetries and compare them to experimental data. We also consider the case where one of the collision partners 
is prolately deformed. Finally, we calculate the excitation function for the total fusion cross section of the system 
i6Q__|_208pJ2 and compare to experiments before to conclude. 
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II. THE TIME-DEPENDENT HARTREE-FOCK APPROACH 

A. Formalism 

Let us first recall some aspects of the TDHF theory. In a non relativistic microscopic approach, the system is 
described by a N-particles state \ip) which is solution of the Schroedinger equation 

ih ^^\^j) = H \^) (1) 

with the microscopic Hamiltonian^ 

i—1 ■i>j^l 



^ We consider only two-body interactions. 



Jdt {ip\H — ih-^\tp) — 0, which is equivalent to Eq. (1), is then solved in the subspace of Slater 



The state \ip) contains all the information on the system, which is more than what we really need for a good description 
of the dynamics. We often need only expectation values of one-body observables, such as the position of the fragments, 
their shapes and particle numbers. These quantities are determined from the one-body density matrix p with elements 
Pij = {tplalailtp) . The expectation value of a one-body observable F = ^- fij alaj is then given by {tp\F\ilj) = Tr (pf) . 
The first step toward the TDHF theory is to restrict the description to one-body observables, and to seek for an 
equation giving the evolution of p. Starting from Eq. (1) and using the Bogolyubov-Born-Green-Kirkwood-Yvon 
(BBGKY) hierarchy [24-26], we can show that the one-body density matrix follows [27] 

ih^^p^[h[p],p]+TT2[vil,2), C(l,2)]. (3) 

where h[p] is the HF single-particle Hamiltonian with matrix elements hij ~ {i\ h[p] \j) — j2_ (^ip\H\ip'^ and C is the 
correlated part of the two-body density matrix. 

Eq. (3) is exact but has two unknown quantities: p and C. The second step toward the TDHF equation is to neglect 
the second term of the right hand side in Eq. (3). This can be done in two alternative ways: 

• The correlation C vanishes if we impose |^) to be an independent particles state at any time. The variational 
principle S 
determinants. 

• The truncation of the BBGKY hierarchy can also be done by neglecting the residual interaction Vres — H — 
'^i^ih[p\{i). This is a mean-field approximation because the Hamiltonian is approximated by a one-body 
operator H ~ ^- hij alaj. In this case, a system described by a Slater determinant at an initial time will be 
an independent particles state at any time. 

We finally get the TDHF equation ih 4^ p — [h[p\ ^ p] where p is now the one-body density matrix of an independent 
particles state. The operator associated to p acts in the Hilbert space of single-particle states. It is written p = 
X]i=i Wi){'^i\ where \ipi) denotes an occupied single-particle state. 

The TDHF theory neglects the pairing correlations which are contained in C. In fact, TDHF describes the evolution 
of occupied single-particle wave functions in the mean field generated by all the particles and assures an exact treatment 
of the Pauli principle during the dynamics. 

B. Practical aspects 

The advantage of TDHF is that it treats static properties and dynamics of nuclei within the same formalism and the 
same interaction. The initial state is obtained through static HF calculations which reproduce well nuclear binding 
energies and deformations. TDHF can be used in two ways for nuclear reactions: 

• A single nucleus is evolved in an external field [8], simulating for instance the Coulomb field of the collision 
partner [28]. 

• The evolution of two nuclei, initially with a zero overlap, is represented by a single Slater determinant [9, 16]. 

The first case is well suited for inelastic scattering, like Coulomb excitation of vibrational and rotational states. 
The second case is used for more violent collisions like fusion reactions. In the latter, the lack of a collision term 
might be a drawback. At low energy, however, fusion is driven by one-body dissipation because the Pauli blocking 
prevents nucleon-nucleon collisions. Fusion occurs by transferring relative motion into internal excitation via one-body 
mechanisms well treated by TDHF. 

Another important advantage of TDHF for near-barrier reaction studies is that it contains all types of couplings 
between the relative motion and internal degrees of freedom whereas in coupled channels calculations one has to 
include them explicitly according to physical intuition. The symmetries corresponding to the internal degrees of 
freedom of interest have to be relaxed in TDHF. However, it gives only classical trajectories for the time-evolution 
and expectation values of one-body observables. In particular, it does not include tunneling of the global wave 
function. 



We use the TDHF3D code built by P. Bonche and coworkers with the SLyAd Skyrme parametrization [17] which is 
a variant of the SLy4 one [29] specifically designed for TDHF calculations^. This code has a plane of symmetry (the 
collision plane). It uses the Skyrme energy functional expressed in Eq. (A. 2) of Ref. [30] where the tensor coupling 
between spin and gradient has been neglected. The step size of the network is 0.8 fm and the step time 0.45 fm/c. 

A TDHF calculation of two colliding nuclei is performed assuming that the two collision partners are initially at a 
distance Dq in their HF ground state. This distance has to be big enough to allow Coulomb excitation in the entrance 
channel (polarization, vibration, rotation...). This initial distance is chosen to be Dq = 44.8 fm. We assume that 
before to reach this distance, the nuclei followed a Rutherford trajectory, which determines their initial velocities Vi 
and V2. The Galilean transformation pi{t = 0) = g"nvir p^p^ ^-imvi-r appjig^ on the HF density matrix of the 
nucleus z (i = 1 or 2) put it into motion with the velocity V; [31]. 

III. FUSION BARRIERS 

Let us consider the simple case of fusion barriers. They are classically defined as the energy threshold above which 
fusion occurs for a head-on collision. Experimentally, the average position of the barrier can be approximated by the 

centroid of the so-called barrier distribution Db (E) — -r— j {(Jfus {E) E) where Ofus is the fusion cross section [32] . In 
the special case of a single barrier for a classical system, the barrier distribution is a Dirac distribution Db ~ 6{E — B), 
whereas a width is generated by tunneling for a quantum system. Several barriers may be generated by the coupling 
between internal degrees of freedom and the relative motion [1]. 

To determine the fusion barrier from TDHF, we consider head-on collisions at various energies. The barrier is then 
located between the highest energy for which there is no fusion and the lowest one for which fusion occurs. The 
average barrier is first studied for spherical nuclei. Then the case of deformed nuclei is considered. 

A. Spherical nuclei 

Fig. 1 shows the density plot for a -'^^0+^°*Pb central collision at a center of mass energy E — 74.44 MeV. After 
a neck formation, the system separates into two fragments. Fig. 2 shows the same reaction at E = 74.45 MeV. We 
see that adding 10 keV is enough to fuse. We deduce the fusion barrier Vb^^^ = 74.445 ± 0.005 MeV. This value is 
in excellent agreement with the experimental one V^^^' ~ 74.5 MeV [33]. It is interesting to note that, if we assume 
frozen HF densities of the collision partners obtained with the same interaction, we get a barrier y^''°^'^" — 76.55 MeV 

at a radius r{^''°'"'"' = 11.73 fm. These values are close to the Bass barrier Vf """ == 77.10 MeV at i?|°"" = 11.42 fm 
[34, 35], but overestimate the experimental fusion barrier. We conclude that the TDHF calculations contain dynamical 
effects that reduce the barrier by ~ 2 MeV as compared to the frozen approximation. It has been suggested that 
transfer may affect the barrier for this system [33]. Indeed, we see in Fig. 2 that the two fragments are linked by a 
neck and form a di-nuclear system during ~ 400 fm/c. It is enough time for the nuclei to transfer nucleons, leading 
to a dynamical evolution of the barrier. 

To get a deeper insight into this transfer preceding fusion, we focus on the case just below the barrier (Fig. 1). Here, 
almost two protons and no neutron, in average, have been transfered from the ^^O to the ^"^Pb. The two-protons 
transfer from the light to the heavy nucleus is then expected to be an important channel at the barrier. This is 
consistent with experimental observations of relatively high Carbon production cross sections, of the same order of 
the Nitrogen ones, in the exit channel of ^^0-t-^°*Pb at the barrier [36, 37]. 




FIG. 1: Densities associated to a ^^O-t-^^^Pb central collision at a center of mass energy E = 74.44 MeV. The surface corresponds 
to an isodensity at half the saturation density. Each plot is separated by 135 fm/c. 



^ The calculations are performed in the laboratory frame and not in the intrinsic frame. It is then necessary to remove the center of mass 
correction in the fitting procedure of the interaction and in the initial HF calculations of the collision partners. The "d" in SLy4d stands 
for dynamics. 




FIG. 2: Same as Fig. 1 at £ = 74.45 MeV. 



Finally, Fig. 3 shows a comparison between experimental fusion barriers and those from TDHF calculations for 
systems with various total masses and mass asymmetries. The lowest barrier is for '*'^Ca+*''Ca and the highest one 
for ^®Ti+^°^Pb. The agreement is of the same order than with the Bass barrier [34, 35]. Remembering that TDHF 
has no adjustable parameter on reaction mechanisms, we conclude that one can use it for fusion barriers prediction 
with confidence. 
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FIG. 3: Theoretical fusion barriers from TDHF calculations (circles) and the Bass barriers (stars) as function of the experimental 
values (from barrier distributions centroid). 



B. Effect of deformation on fusion 

We now consider collisions of a spherical nucleus on a deformed one. In such a case, the barrier depends on the 
orientation of the deformed nucleus at the touching point, leading to a wider barrier distribution than the single 
barrier case [1]. 



Heavy deformed + light spherical nuclei 



Fig. 4 shows barrier distributions for collisions with a light spherical projectile (^^0) on heavy prolately deformed 
targets ^^^Sm (left) and ^ss^j (right) [38, 39]. Their width are ^7—10 MeV and cannot be explained by tunneling 
alone. The latter adds a width of only 2-3 MeV [32] . Such barrier distributions are usually well reproduced in the 
framework of coupled channel calculations [1] . In addition, microscopic theories like TDHF can help to understand 
the physics process generating these couplings. 

The barriers predicted by the TDHF calculations for these systems are also shown (arrows) on Fig. 4 for two extreme 
configurations of central collisions where the collision and deformation axis are either parallel or perpendicular. In 
the parallel configuration, the Coulomb repulsion is smaller and so is the resulting barrier. All the intermediate 
orientations at the touching point give fusion barriers between the parallel and perpendicular configurations ones. 
In addition to a good reproduction of the centroid, the TDHF calculations also reproduce the width of the barrier 
distributions generated by a static deformation of the target without any adjustment of parameters. 
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FIG. 4: Experimental barrier distributions for ^''O+^^^Sm (left) and ^''0+^^**U (right). The arrows indicate the barriers 
obtained from TDHF calculations for two extreme configurations where the collision and deformation axis are parallel (lowest 
barriers) and perpendicular (highest barriers). 

The shape of these distributions is due to a prolate deformation of the target: the distribution is more peaked 
at higher energies. The configuration where the coUision and deformation axis are parallel corresponds to only 
one possible orientation, whereas the perpendicular configuration can be reached by any orientation for which the 
deformation axis is contained in the plane perpendicular to the collision axis. If one assumes an isotropic distribution 
of the orientations at the touching point, which is a reasonable approximation when the spherical nucleus is light and 
the deformed one heavy [28], then the perpendicular configuration is more probable than the parallel one and the 
barrier distribution is peaked at high energies. One gets the opposite with an oblate nucleus instead of a prolate one. 

2. light deformed + heavy spherical nuclei 

The case of a light deformed projectile on a heavy spherical target has been investigated theoretically in Ref. [28] 
both within the TDHF and the coupled channel frameworks with the code CCFULL [40]. For such systems, the 
barrier distribution gets affected due to reorientation of the deformed nucleus in the Coulomb field of the target. This 
breaks the isotropy of the orientation axis distribution and results to a fusion hindrance at low energies. Experimental 
evidences of this effect have been reported recently [41]. The reorientation is proportional to Asphericai / Atotai and 
then does not affect the systems studied in Fig. 4"^. 

IV. EXCITATION FUNCTION OF ^^^O+^'^PB 

We now focus on fusion cross sections given by 



t2 °° 
^ 1=0 



(4) 



where /i is the reduced mass of the system and Pfus{E, I) is the fusion probability at a center of mass energy E and 
an angular momentum \/l(l + l)h. The restriction to an independent particles state, as in TDHF, leads to Pfus = 1 
for I < Imax and for / > Imax- We then get the so-called "quantum sharp cut-off formula" [42] 



<7fus{E) 



2fiE 



{.'"max ~F ij • 



(5) 



To avoid discontinuities introduced by the cut-off and the integer nature of Imax, we approximate (Imax + l)ft by its 
classical equivalent Cc- The latter is the threshold of the classical angular momentum C = ^/2jlEb, b being the impact 



^ Though it is a Coulomb effect, the reorientation does not depend on the charges of the nuclei but on their masses as we can see in 
Eq. (5) of Ref. [28]. 



parameter, below which fusion occurs [35]. This approximation is justified by the fact that both {Imax + 1)^ and C'^/h'^ 
are greater than lmax(Jmax + 1) and smaller than (Imax + ^){lmax + 2). We finally get the classical expression for the 
fusion cross section afus{E) ~ TT£'^/2fj,E. 

Figure 5 shows the excitation function obtained for the ^^O+^^^Pb system in comparison to experimental data 
[33]. There is a good agreement above the barrier, though the fusion cross sections are overestimated by about 16 %. 
However, the fusion cross section vanishes below the barrier, following a classical behavior. This drawback of TDHF 
is well known and is due to the restriction to a single independent particles state. Indeed, to get a fusion probability 
between and 1, we need at least two Slater determinants: one describing the two well separated fragments after the 
collision when fusion does not occur and one describing the fused system. It is necessary to go beyond TDHF to treat 
a sum of Slater determinants and then to describe sub-barrier fusion due to tunneling effects. 
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FIG. 5: Experimental fusion excitation function (stars) for 
fusion cross sections obtained from TDHF calculations. 



3Q^208pj^ The lines denote the upper and lower limits for the 



V. CONCLUSIONS AND PERSPECTIVES 

We presented some applications of the Time-Dependent Hartrec-Fock theory to nuclear fusion. The only phenomeno- 
logical input is the set of parameters of the SLy4d force which have not been adjusted on any reaction mechanism 
like cross sections or fusion barriers for instance. Despite this, the agreement between the TDHF calculations and 
the experimental fusion barriers is excellent for a wide range of projectiles and targets. The width of fusion barrier 
distributions generated by a static deformation of the target is also well reproduced. 

The TDHF calculations overestimate the fusion cross sections for the system ^®0+^"^Pb above the barrier by 
about 16%. Below the barrier, the fusion cross section vanishes in TDHF calculations. The sub-barrier fusion due to 
quantum tunneling of the many-body wave function is not present in TDHF. This is due to the restriction to a single 
independent particles state. 

Though the fusion cross section has a classical behavior, the quantum nature of the single-particle wave functions 
is well treated. Thus, it would be interesting to study the transfer of nucleons from one nucleus to the other within 
TDHF''. An illustration of such transfer is seen in the TDHF calculation of the ^^0-|-^°*Pb just below the barrier. 
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